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Abstract 

Experiments on quasistatic crack propagation in gelatin hydrogels reveal a new branching insta- 
bility triggered by wetting the tip opening with a drop of aqueous solvent less viscous than the bulk 
one. We show that the emergence of unstable branches results from a balance between the rate 
of secondary crack growth and the rate of advection away from a non-linear elastic region of size 
Q/E where Q is the fracture energy and E the small strain Young modulus. We build a minimal, 
predictive model that combines mechanical characteristics of this mesoscopic region and physical 
features of the process zone. It accounts for the details of the stability diagram and lends support 
to the idea that non-linear elasticity plays a critical role in crack front instabilities. 

PACS numbers: 46.50.+a,62.20.mt 
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Introduction 

Recent developments in tissue engineering |l| have raised polymer hydrogels to the status 
of genuine, structural materials, suitable for load bearing applications, such as scaffolds for 
in vivo tissue regeneration. Disregarding the biological issue, an essential step towards the 
rational design of such hydrogel-based systems is to understand how these "soft" solids break. 
This means not only relating their micro-structural features to their mechanical strength, 
but also describing the way failure will eventually proceeds, i.e. how damage will pervade a 
sample set in a given mechanical as well as physico-chemical environment. 

Hydrogels are network-based, solvent-swollen materials. The covalent (resp. non cova- 
lent) nature of the cross-linking bonds plays an essential role in chemical (resp. physical) 
gel strength. Chemical gels require C-C backbone scission to break and exhibit quasi rate- 
independent energy dissipation during crack propagation. Hence, fracture usually proceeds 
dynamically, i.e. accelerates until reaching values in the vicinity of the Rayleigh wave ve- 
locity of the material. Chemical gels are therefore termed brittle and as such have received 
much attention recently for their ability to mimics the dynamic fracture phenomenology 
of hard brittle solids, albeit at much lower velocities This is in contradistinction with 
physical gels in which weaker crosslinks act as mechanical fuses preventing chain scission. 
Disrupting large crosslink structures and/or pulling chains out of the gel matrix ^\ are 
usually strongly dissipative, rate-dependent breaking processes which are responsible for the 
amazing resistance to fracture of these materials, mostly consisting of solvent with only a 
small fraction of polymer. Steady, quasistatic crack propagation is readily achieved with 
physical gels. This makes them suitable for detailed investigation aimed at unraveling the 
physical, dissipative mechanisms at work. 

Notwithstanding these fundamental differences, both physical and chemical gels have 
in common their softness i.e. their large elastic compliance. Consequently, the fracture 
process involves finite crack tip opening displacements and large deformations over a wide 
region ahead of the tip. This situation is markedly at odds with the requirement of linear 
elastic fracture mechanics (LEFM) that small strain, linear elasticity prevails everywhere 
apart from a near tip "process zone" where dissipative bond-breaking occurs. Owing to the 
powerful ability of LEFM to predict the onset of failure with minimal material-dependent 
input js], though, the role played by the near crack-tip, non-linear elastic (NLE) zone where 
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LEFM breaks down, has been long overlooked. This, nevertheless, cannot be ignored any 
more when dealing with two important issues: 

Elastic crack blunting — In soft elastic solids which can sustain bond-breaking stresses 
much larger than their small-strain Young modulus E, crack tips tend to blunt, i.e. to 
develop large radii due to purely elastic deformations [6], thereby mitigating remote stress 
intensification. A strict LEFM approach would lead to the paradoxical conclusion that 
stresses, which are predicted to plateau at values of order E, remain too low to initiate 
rupture Gj. Strain-hardening, a NLE characteristics shared by strongly stretched polymer 
networks, has been proposed as a likely way of restoring stress concentration in the near tip 
region. 

Crack front instabilities — Smooth crack propagation, either quasistatic or dynamic, is 
rather the exception than the rule in fracture of soft elastic solids. Branching j^, splitting 

m 

or oscillating cracks are commonly reported phenomena which LEFM fails to 

explain Again, the existence of the NLE zone, introducing a new length scale in the 
fracture problem, has been invoked as the missing ingredient for predicting the onset of a 



t is worth noting that crack tip 



13| as possible causes of the high 



front instability and its characteristic features [UJ, Il2|. 
blunting and splitting have already been evoked by Gent 1 
tear strength of visco-elastic elastomers (unswollen chemical gels). He listed them amongst 
several unresolved issues in rubber fracture. Indeed, despite recent advances in numerical 
and theoretical description of the NLE zone, elastic blunting and front instabilities remain 
widely open problems which probably transcend the case of hydrogels. 

These materials, which feature no noticeable linear visco-elasticity over a wide frequency 
range, are not expected to exhibit bulk energy dissipation during crack propagation but ra- 
ther localized dissipation in a near tip process zone of extension d. In the case of gelatin gels, 
it has been argued, based on experiments [ij] that c? ~ 100 nm. The size of the NLE zone 
scales with the natural length in fracture problems C = Q/E where Q is the energy release 
rate (free energy released by the advance of a unit area of crack) which, for a quasistatic 
crack, identifies itself with the dissipated fracture energy. In the case of gels, which are 
"tough" solids with relatively large Q and low E, C ranges typically between 100 fim and a 
few mm (as compared, e.g. to a few A for brittle silica glass). Thus, fracture in hydrogels 
exhibit a clear hierarchy of relevant length scales between the microscopic d, the mesoscopic 
C and the macroscopic system size h : 
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This configuration, which extends the so called "small scale yielding" hypothesis la] to the 
NLE case, plays a central role in the following. 

This article aims at presenting experimental evidence of a previously undescribed branch- 
ing instability in a physical gelatin gel. The clear separation of length scales enables us to 
unravel NLE effects from dissipative mechanisms. The instability is therefore amenable to 
a physical interpretation which accounts — though schematically — for the elastic blunting 
of the quasistatic main crack on the scale C More precisely, the unstable onset of side 
branches, which we trigger by modifying dissipation locally, i.e. within the process zone 
itself, is proposed to result from the competition between the growth of secondary cracks 
and their advection in the displacement field of the main crack, out of the near tip NLE zone 
where opening stresses remain significant. 

The article is organized as follows: in section [U the physical mechanisms of gelatin 
gel fracture are briefly outlined, with emphasis put on the role of solvent viscosity. The 
branching instability experiment is described in section [TTl and the marginal nature of the 
branching onset is evidenced. This makes it possible in section IIIII to propose a model 
featuring a single free parameter, lumping together the NLE properties of the blunted crack. 
Once this parameter is determined, the predictability of the model is successfully tested. 
Possible generalization of the model to other soft solids is discussed in part | 



I. A REMINDER ON GELATIN GELS 



A. Structure 



Gelatin, a biopolymer made of denatured collagen, can dissolve in aqueous solvents (e.g. 
mixtures of water and glycerol) above the melting temperature (T^ ^ 40°C). Upon cooling 
below gelatin chains partially revert to the triple helix conformation of the native colla- 
gen, resulting in a network of rigid rods interspeded by random coils [l^]. The mesh-size ^ 
can be evaluated from the Young elastic modulus E assuming an entropic origin for the gel 
elasticity: ^ ~ (ksT/Ey^^. It is typically on the order of 10 nm for Young moduli in the 10 
kPa range. 
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The cross-linking triple helices are stabilized by weak physical H-bonds. Accordingly, the 
gel is thermoreversible. Moreover, it ages as revealed by the logarithmic increase of E with 



time, and exhibits slow stress relaxation typical of soft glassy materials jl8l]. However, on 
short times relevant to most fracture experiments, the polymer network can sustain shear. 
The relative motion of the solvent with respect to the elastic network is a collective, diffusive 
process with a coefficient Dcoii ~ -E'^^/''7buik with ?7buik the solvent viscosity. Typically Dcoii — 
10^^^ m^.s~^, an order of magnitude which makes solvent draining a very slow process 
on macroscopic length scales so that, for all practical purposes, the gel samples can be 
considered quasi-incompressible. 

Non-linear elasticity of gelatin gels exhibits strain hardening at moderate strains. As 
shown on figure [T], it does not obey the classical neo-hookean constitutive law (cor- 
responding to a strain energy density functional Wnh{^x, ^y, ^z) = EJi/6 with Ji = 
A^. + + — 3 an invariant function of the principal stretch ratios, related by the in- 
compressibility condition A^jAyA^ = 1). The data obtained in uniaxial compression are, 
rather, consistent with an empirical expression used in numerical fracture studies as a sim- 
ple model for severe strain hardening pj!]: Wsn = EJm{exp{Ji/Jm) — l)/6 with Jm — 2.3 
suggesting that for elongations of order + 3 ~ 230% the chains are already stretched 
significantly taut. 

B. Fracture 

Experimentally, the rate-dependent fracture energy Q{V) for a crack propagating steadily 
at a markedly subsonic (or "quasistatic" ) velocity V reads : 

g = Qo + rr/buikl^ (1) 

with r^buik the solvent viscosity, — 1 J.m~^ and F ~ 10®. This linear relationship has been 
interpreted ^ as resulting from the fact that in this physical gel, fracture does not proceed 
by chain scission but rather via unzipping of the cross-links and subsequent pulling out of 
the overall chains to the expense of the viscous drag against the solvent. Qq stems from 
the plastic work for unzipping and, possibly, from the "dehydration" cost for exposing the 
polymer chains to air. The rate-dependent term accounts for the viscous losses against the 
solvent, r is predicted to scale as the squared ratio of the chain contour length (A ~ 1 fim) 
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FIG. 1: Non-linear elastic response of a cylinder of gelatin gel (5 wt.% in 60%-40% glycerol-water 
solvent) in uniaxial compression between well lubricated plates, a is the true (Cauchy) stress and 
Xx is the compression ratio. The solid curve is a fit using a strain energy density Wsk (see text) 
with Jm = 2.3. The initial slope (A^. < 1) yields the Young modulus E = 12.1 kPa. The dashed 



line corresponds to the neo-hookean elastic solid with the same small strain modulus. The sma. 



hysteresis between loading and unloading curves is probably indicative of slow stress relaxation 
during the 6 s duration of the cycle. 

to the mesh size of the network ~ 10 nm). 

Note that the stress level at the crack tip, given by ^/A is always several order of mag- 
nitude larger than the Young modulus. This is precisely the case where elastic blunting of 
the crack is predicted, hence where a marked NLE zone ahead of the crack tip is expected 

0. 

The sensitivity of the fracture energy to solvent viscosity makes it possible to estimate 
the extension d of the process zone where energy is dissipated [4]- This is achieved by 
wetting the crack tip opening with a drop of solvent of viscosity ?7drop < Vhuik- Both drop 
and bulk solvents differ only by their fractions of glycerol in water which tend to equilibrate 
via molecular diffusion, with a coefficient -Dgiy. Experimentally, the drop acts as a reservoir. 
The glycerol content of the gel ahead of the moving tip adjusts from the drop concentration 
to the bulk one over the diffusive skin depth of order Dgiy/V. At low enough velocities, the 
effective viscosity 



Vcs = r 



dV 



(2) 



which accounts for the advance of "rinsing" within the process zone is that of the drop 
liquid. As V increases above Vdifr — Dgiy/d, ries{V) gradually grows up to r^buik- Measuring 



Vdiflf yields li ^ 100 nm. 



This estimate was obtained originally [4] with a relatively small viscosity contrast, 
TponYk/Vdrop — 3. Unexpectedly, on attempting to reproduce the experiment with a larger 
contrast, a spectacular branching instability was observed. It is the aim of this paper to 
describe and interpret this phenomenon. 

II. EXPERIMENTAL 

A. Material and methods 

Gel samples are composed of 5 wt.% of gelatin (Sigma, "300 Bloom" grade) in aqueous 
solvents containing from = to 70 wt.% of glycerol in water. The pregel solution is 
prepared by letting gelatin dissolve in its solvent at 85° C under gentle stirring. It is poured 
into rectangular molds (length L = 300 mm, width h = 30 mm, and thickness e = 10 
mm) and "set" at 5°C for 10 hours. Before performing any experiment, the samples are 
left to equilibrate at room temperature (T ~ 20°C) for 2 hours. Aging during the duration 
of an experimental run (about 2 min) has been checked to be negligible. The samples are 
characterized by their small strain Young moduli E (between 9 and 15 kPa, increasing with 
glycerol content) and their solvent viscosity r^buik (from 1 cP for pure water to 22 cP at 



FIG. 2: Schematic experimental setup. The gel sample is stretched between two rigid, parallel grips 
so that the energy release rate is imposed. The steady-state crack front velocity V is measured. 

Wetted crack experiments are performed with a setup fully described elsewhere [l^ . Gel 
plates are stretched along their width by imposing a displacement Ah to the rigid grips 
(see Figl2]). Due to the large aspect ratio L/h = 10, uniform stretching is realized in a 
significant portion of the sample, thereby storing a prescribed amount of elastic energy Q to 
be ultimately released per unit area of a straight crack propagating along the length. Q ( A/i) 



= 70%). 
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is evaluated as W/(eL) with W the measured work required to stretch an un- notched sample 
wAh. Note that, although the expression for Q is strictly valid for infinitely long samples 



151], it is expected to provide a reasonable approximation to within about h/L = 10 ^. 



The loading rate for the calibration test is chosen high enough for stress relaxation to be 
negligible during the load (see eg. Fig[T]). Also the phenomenon reported in the present 
study occurs on time scales short enough for internal dissipation associated with remote 
stress relaxation to be negligible in comparison with the energy released by the fracture 
itself. 

Crack propagation is initiated by cutting a notch at one end of the plate. Away from 
the sample edges, when the crack velocity reaches a steady value V of order a few mm.s~^, 
a drop (~ 200 yuL) of a water /glycerol mixture of viscosity //drop is quickly injected into the 
tip opening where it remains trapped by capillarity and gravity (the crack travels downward 
vertically). The subsequent dynamics of the crack tip is monitored at 10 frames per second 
by a video camera mounted on a traveling stage so as to keep the crack front within the 
field of view. 



B. Characteristics of the branching instability 

Figure [3] shows the first stage of the instability triggered by pure water (r^dmp = 1 cP) 
wetting the tip of a gel with cj) = 60% (r/buik = 11 cP). Within a few seconds, the gel ahead 
of the crack is pervaded by a 3D damaged zone made of microcracks. The zone coarsens 
until the average front becomes flat. At that point, the fracture process is almost inhibited 



(for the subsequent slow return to a straight, dry crack, see movie [20|). 

This spectacular, bursting response to a rather modest environmental change is in marked 
contrast to our previous report [14] on the wetting by pure water of a = 30% gel (r/buik = 1.8 
cP). Then, while being markedly accelerated the crack remained straight. Crack acceleration 
was ascribed to two solvent effects: firstly, wetting the tip prevents the extracted chains from 
being exposed to air, hence a significant lowering of the threshold energy Qq; secondly, the 
water drop induces an osmotic imbalance, hence a glycerol depletion ahead of the tip, as 
described in II B[ Both effects contribute to cutting down the fracture energy. Since the 
energy release rate is imposed, it results in the observed speeding up of the crack. FigO] 
shows that, beyond this effect, local solvent dilution can destabilize a crack tip as well. 
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FIG. 3: Snapshots (At = 0.3 s) extracted from a movie [2d] showing the coarsening of a microc- 
racked zone ahead of a crack tip wetted at t = by a drop of pure water in a = 60% gel. Crack 
tip opening is 12 mm {Q = 65 J.m~^). The graph displays the front velocity before and after tip 
wetting. 



depending on the bulk viscosity. 
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FIG. 4: Stability diagram at constant crack tip opening (12 mm) for cracks in gels with bulk 
water /glycerol solvents of viscosities ?7buik wetted by drops of viscosities r/drop- The shaded zone is 
out of experimental reach. Near critical data are labeled by the values of the compound parameter 
(%/%uik)-(^/^dry) with V and Vdry defined on Figl21 The dashed line is the predicted locus of 
critical cracks (see text, section UllPp . 

In order to clarify this point we have mapped the stability of the crack onto the 
(^uik, ''^drop) plane for a prescribed 12 mm crack tip opening. As displayed on FigJH branch- 
ing occurs for large enough viscosity contrasts. Since the tip solvent cannot be made more 
inviscid than pure water, there is a minimal bulk viscosity, corresponding to (p = 50% below 
which the crack remains stable whatever the tip environment. Figj5] shows a snapshot of a 
crack in this marginal gel (0 = 50%), 4.8 seconds after being wetted by a drop of pure water. 
Though no thick damaged zone develops ahead of the tip, the crack path exhibits distinct 
undulations and aborted side branches. Moreover, before these secondary cracks stopped, 
they started themselves to branch, which we interpret as an interrupted cascading process 
that, if unimpeded, would have led to the microcracked zone of FiglH] 

Now, for a gel with a bulk viscosity above the marginal one, increasing the drop viscosity 
from that of pure water changes the state of crack propagation from the unstable one shown 
on Figl3]to a stable one via a regime similar to that described above (see movies [2o|). It is 
therefore legitimate to term "critical" this regime of crack emission which marks systemat- 
ically the frontier between stable and unstable fracture. 
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FIG. 5: Wet crack in a critical state exhibiting aborted side branches. 



III. A SCHEMATIC MODEL 

Disregarding the complex, fully developed structure of FigjSl we now focus on the critical 
regime and propose a criterion for the onset of instability that combines the microscopic 
model of fracture energy specific for physical hydrogels such as gelatin and the more generic 
concept of a mesoscopic near-tip, non-linear elastic (NLE) zone. As a starting point we 
note that secondary cracks are emitted in the close vicinity of the main tip. While they 
grow, they are progressively advected away and gradually close on approaching the straight, 
traction free crack edges. Thus, there is a competition between growth and advection rates. 
We propose that only those cracks which manage to grow significantly before being advected 
out of the near-tip zone, where opening stresses are significant, can in turn branch, hence 



trigger a cascading process. We therefore describe branching as a convective instability 2^ ■ 



A. Near-tip stress field of a blunted crack 

Branching requires nucleation and growth of secondary cracks which eventually compete 
with the initial main crack to release elastic energy. We assume in the following that nucle- 
ation sites, or flaws, are provided by the inhomogeneities — of either toughness or stiffness 
— of the randomly crosslinked polymer network. 

As a consequence of elastic blunting, associated to strain-hardening of the gel, the region 
where opening stresses are high is along the crack face rather than directly ahead of the 
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crack tip, a qualitative feature which is clearly revealed by finite element simulations 12|. 
This suggests that branching probably originates from nucleation and growth of edge cracks 
in this quasi-uniaxially stressed zone. Its extension is provided, in a first approximation, 
by the radius of curvature of the near tip parabolic crack opening. This is given, within a 
multiplicative constant of order unity, by Q / E. This classical result of LEFM ^] remains 
valid in the NLE case, at least for a neo-hookean solid (see e.g. |l2]), and is certainly 
dimensionally correct for more strongly strain hardening materials. 

We now aim at writing a growth equation for an edge crack in the NLE stress field. This 
is a desperately complex mechanical problem. We therefore make one further step and, 
based on the previously described structure of the near-tip stress field 12|, we schematize 
the blunted elastic zone as a strip of size ~ G/E, uniformally stretched by a stress ctq 
perpendicular to the fracture plane (see Figl6]). The material is assumed to have strain- 
hardened up to a state characterized by a uniform small strain modulus E^.^ > E. 

Progressive building-up of strain hardening on approaching the blunted crack tip is now 
replaced by a stepwise jump in small strain modulus at a distance Q/E. The unknown stress 
(Jo is determined according to the compatibility requirement that LEFM should hold in the 
outer region. There, the stress intensity factor [5] Kf ~ VGE, hence the stress at a distance 
Q/E is (Jo ~ Ki/^g/E ~ E. Note that this rough schematization of the blunted elastic 
zone is akin to that proposed in Jg] on the basis of slightly different arguments. 



B. Growth mechanism for a secondary crack 

Since we are interested in the onset of unstable branches, we assume the length ^ of the 
secondary crack to remain small as compared to the extension Q / E oi the blunted crack zone, 
so that it can be treated as an edge crack which does not interfere significantly with the stress 
field of the main crack. For the sake of simplicity, we estimate the corresponding energy 
release rate ^branch according to LEFM [sl in a material with an effective Young modulus 
E'eff. Within a factor of order unity: ^branch = (7r(Jo/-£'ofr)^ = PEE, where we have introduced 
a dimensionless factor (3 ~ E/E^s which is indicative of the level of strain-hardening. The 
clear separation of lengthscales, d <^ Q/E, makes it legitimate to assume that the edge crack 
dynamics on scale Q/E is ruled by the same physical mechanism as the main crack. The 
growth equation is thus obtained by equating ^branch with the fracture energy Q given by 



12 



([T]), replacing V by d£/dt: 



d£ 

pEe = go + r?7branch3T 

at 



(3) 



where r/branch is the effective solvent viscosity "felt" by the chains pulled out of the process 
zone ahead of the secondary crack. 

As discussed in the following section, a possible rate-dependence of r/branch can be ruled 
out. Solving eq. (jl]) for i(t) is therefore straightforward : £(t) = ic + (4 — ^o) exp(i(:/r). Ac- 
cordingly, any supercritical seed crack of initial size £o > = Qo/ (3E will grow exponentially 
over a characteristic time 

r ^branch 



r 



(3 E 



(4) 




i i i i 



■Q/E 



FIG. 6: Schematic representation of a secondary crack growth and advection in the NLE zone of 
a blunted main crack. 

If the secondary crack is left to grow in the constant stress field it will eventually, at times 
on the order of r, compete with the main crack for releasing energy, and give rise, in turn, 
to side branches. However, while the main crack grows, the root of the secondary crack is 
advected away from the main tip, eventually leaving the NLE zone where sufficient stress 
exists to permit branch growth. The stability criterion discussed in the following is based 
on the competition between these two effects. 
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C. A criterion for stable, wetted cracks 



The life time of a secondary crack is given by the duration of its travel along the blunted 
tip zone as it is advected away from the apex of the main crack. An upper value for this 
duration is the time required for crossing a region of size Q/E at V, the velocity of the main 
crack tip in its wetted state : 

T^adv = (5) 

As long as r > Tadv, all secondary cracks will stay in their embryonic states and never 
compete with the main crack which will therefore remain stable. Since, according to ([2]), 
is given implicitly for a prescribed Q by 

g = go + rr],s{v)v (6) 

with Qq being usually negligible (see FigJTj), the stability criterion reads: rj^^ranch/Vcsiy) > P- 
At this stage, it is crucial to realize that the dynamics of the secondary crack is ruled 
by the effective viscosity r/branch which depends on glycerol diffusion, i.e. on the velocity of 
the branch tip with respect to the liquid drop. Since the latter is driven by capillarity and 
gravity along with the main crack, the relevant velocity is that at which the secondary tip 
cuts into the primary crack lips, i.e. di/dt. We have therefore ^^branch = Vesidi/dt). Since 
di/dt remains much smaller than the velocity V of the main tip during the maturation phase 
(secondary crack growth times t < r), one can expect that r]es{di/dt) < rics{V). In fact, it 
can be shown that r/branch = "/^drop (see appendix) so that, finally, the criterion for a crack to 
remain stable at a quasistatic velocity V reads simply: 

>P (7) 



^eff(V^) 

D. Experimental tests of the model 

The above stability criterion ([7]), is expressed in terms of the effective viscosity func- 
tional, defined by ([2]), and parameter (3. The former accounts for the local modification of 
the microscopic process zone via solute diffusion, the latter is aimed at catching the NLE 
modification of the gel at the mesoscopic scale of the blunted crack opening. 

Let us note in the first place that, in the destabilizing configuration where the drop is less 
viscid than the bulk solvent, the inequality V: rjdrop/VesiV) < 1 holds for any combinations 
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of V, ?7drop and r^buik, including those corresponding to stable cracks. Hence, according to 
([7]), /3 < 1. So, in spite of the numerous multiplicative constants of order unity hidden in /3, 
this parameter therefore retains its physical flavour as a qualitative indicator of the level of 
strain hardening : f3 ~ E/Ecs < 1. 

Furthermore, since rj^s is an increasing function of the tip velocity, ([7j) predicts that 
cracks are stable at low enough velocities. This is a counter intuitive result, owing to the 
proposed convective nature of the instability, since for a slow main crack, side branches are 
only slowly advected away from the tip zone, hence have a priori more time to grow. The 
fact that this direct effect of V disappears in ([71) can be traced back to the V^-dependence 
of the NLE length scale Q{V)/E. A slow crack offers a shorter NLE zone for branches to 
grow. The existence of a lower critical velocity for a given pair r^jrop < ''7buik is therefore a 
strong test of the validity of our schematic model. We have checked it by increasing stepwise 
the crack opening, hence Q, while keeping the tip wet by continuous solvent dripping. FigIT] 
shows the results for a pure water drop wetting a = 60% gel. Indeed, the crack remains 
stable as long as V < Vc with Vc — 4.7 mm.s~^. 

80 
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FIG. 7: Stability diagram in the fracture energy vs. crack-tip velocity plane. Unwetted (□) cracks; 
stable (O) and unstable (•) cracks wetted by pure water (for unstable cracks, V is the front velocity 
measured just after wetting, see Figl2|). The curve is a guide for the eyes. Slopes of the thin, thick 
and dashed lines are in the ratio r/tip : rjcs{Vc) : ?/buik- Vdm and Vc are defined in the text. 

This enables us to evaluate (3 quantitatively. Following the above discussion, for V < Vdig 
the tip is fully rinsed and rj^s — ''7drop, hence Q ~ "/^dropr^- Fitting the low velocity data 
accordingly yields F ~ 1.6 x 10^, a value compatible with that of non- wetted cracks in a 
companion sample of the same gel (see Fig. [7]). At the critical point, we measure ?7efr(K;) = 
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{Q{Vc) - Go)/O^Vc) from which we deduce that (3 ~ 0.42. 

Finally, we return to the stability diagram of FigJH Though chosen for the sake of its 
experimental simplicity, the constant crack tip opening prescription is theoretically awk- 
ward since it corresponds to a non trivial section of the critical surface in the 3D control 
parameter space (//bulk, "/^drop, It is nevertheless possible to make use of this diagram 
for an independent check of the model by recasting criterion ([7]) in terms of Vjry (resp. 
V^), the crack velocities before (resp. just after) tip wetting (see Figj3]). Neglecting Qq, 
the constant opening imposes Q ~ r/buikrVdry — Ve{i(y)^V, hence the stability condition : 
('7drop/''7buik)-(V^/Vdry) > P- We havc computed the left hand side compound parameter 
for near critical data on both sides of the bifurcation (see FigHj). Remarkably, the figures 
systematically bracket the value f3 = 0.42 determined from the critical velocity. 



IV. CONCLUDING REMARKS 



We have identified and analyzed a new branching instability occurring in the tip vicinity 
of a quasistatic crack in gelatin hydrogels. We have proposed what we think is the first 
predictive model of a convective crack branching instability in a disordered elastic material. 
The key ingredient is the existence of a mesoscopic zone ahead of the crack where residual 
opening stresses are large enough for promoting side branching. The finite size of this near 
tip zone, which scales with Q/E, results in a finite advection time of incipient branches 
away from the main crack. Whether catastrophic crack growth occurs during this time 
chiefly depends on the rate-dependence of the fracture energy. In the case of gelatin, it can 
be conveniently tailored by a local "osmotic" control of the process zone. 

Beyond the deciphering of a puzzling phenomenon which might, at first sight, appear as 
an idiosyncrasy of gelatin gels, we think that our analysis lends support to the emerging 
idea that the non-linear elastic field which bridges between the process zone and the generic 
linear elastic region, plays a crucial role in crack dynamics and, more specifically, in crack 
front instabilities [22|. Though branching instabilities have been mostly documented for 
fast cracks, our study indicates that they could prevail as well for slow ones over wide do- 
mains of their mechanical control parameters (crack tip velocity, tip environment, structural 
disorder,. . . ). 

To go one step further along this line, let us analyze a case where quasistatic crack 
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propagation would be ruled by a fracture energy that depends as a sublinear power-law on 
the crack velocity, say: 



g = go[l + [V/VoT] with a < 1 

Such a functional form is typical of elastomers [l^ and has been reported for gels made of 
solvated triblock copolymers as well sl- In the case of elastomers, the material-dependent 
reference velocity Vq depends on temperature, presumably via the visco-elastic spectrum 
of the material. We assume that the generalized small-scale yielding approximation {d ^ 
g/E <^ L) holds here. 

The growth equation g{i) = jSEi for an edge crack nucleated within the primary blunted 
tip region can be solved analytically. It is straightforward to establish that the degree of 
supercriticality of the crack, introduced for the sake of simplicity as e(t) = — 1 with 

= go/PE, reads: 

e(t) = ^^-i^withr' = ^and7 = a-^-l>0 

Crack length divergence at finite time r' signals an instability. It is important to note, 
however, that the time r' for catastrophic failure depends on eo = e(0), the initial degree of 
supercriticality of the unstable crack seed. 

Due to advection of the material in the tip region, the main crack will however remain 
stable against side branching provided that Vt' > g(y)/E, i.e. 

""/^^ > P7e2 (8) 



1 + {v/Voy 

Here, in contradistinction with the a = 1 case, a secondary crack will close back before 
exploding if the main one is driven at a high enough velocity (the rhs of (IHl) is an increasing 
function of V/Vq). Yet, this upper critical velocity depends on the strength eo of the nu- 
cleation seed. The larger the flaw, the more destabilizing it is for the main crack. Such an 
instability would possibly manifest itself as follows: in the absence of any external triggering, 
a crack running steadily at a velocity V would eventually become unstable, provided it meets 
a flaw of sufficient strength. Since, according to (jH]), the size of the critical seed increases 
with V, the higher the velocity, the lower the probability of meeting such a nucleation site 
and the longer, statistically, the crack would remain stable. 
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This example shows that, even in its simplest version, our model predicts a wealth of 
qualitative behaviors. It prompts us to have a fresh look at the complex, intermittent 
dynamics reported in the case of triblock copolymer gels js], and still unexplained. It also 
points to the important role played by the structural disorder on branching and suggests to 
study controlled, inhomogeneous materials such as filled elastomers or bicomponent, phase- 
separated polymer hydrogels. 



Appendix 

In this appendix, we show that the velocity of a branch is, in its maturing phase, small 
enough so that its process zone is fully rinsed by the drop i.e. that r/branch = ''7drop holds in 
eq. H 

Consider a supercritical crack seed of initial length i{0) = (1 + eo)£c- For t < t, the 
crack velocity remains close to its initial value vq = eo^c/T = eo^o/(^branchr)- Since rebranch > 
^?drop > ^/water, and — 1 J.m~^, t>o < cq X 1 mm.s~^. We argue in the following that, in 
our experiments, eo remains small enough so that Vq < Vdig. 

Our model assumes that "flaws" preexist is the sample and act as edge crack seeds when 
reached by the blunted tip zone of the main crack. As already mentioned, flaws are ascribed 
to frozen fluctuations of the polymer network structure e.g. via its crosslink density. 
Such inhomogeneities remains after coarse-graining up to the scale d of the process zone and 
result in a spatially modulated gel strength about an average value Qq, as evidenced by the 
roughness of fracture surfaces in gelatin gels Q]- An estimate of the maximum amplitude 
of Qo variations, AQq/Qq ~ 14%, has been obtained by analyzing the pinning of very slow 
crack fronts induced by "tough" sp^ts [3]. 



According to Griffith's theory 



3, 



crack initiation results from a balance between the 
fracture energy cost Qo and the elastic energy released in a deformed region scaling with 
i. That is a nucleation process, with an activation energy corresponding to a critical crack 
length ic- With — 1 J.m~^, £c = Qq/[3E ~ 100 ;um. The activation barrier £^act ~ ^o4e, 
with e the lateral extension of the crack, is always orders of magnitude larger than ksT 
whence thermal activation cannot be responsible for the distribution of initial crack lengths. 
We propose here that it rather reflects the statistical distribution of the flaw strengths Qq, 
at the scale of d and therefore identify the supercriticality index eo with A^o/^o so that 
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vo<l mm.s-^ With -Dgiy(0 = 60%) = 2 x 10"^° m^.s"^ and d ~ 100 nm, Vdis ^ 2 mm.s-^ 
Since vq < Vdis, one may therefore safely consider that jybranch — Ves{vo) — Jydrop- 
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